> 

o 



October 1994 Wash. U. HEP/94-63 

hep-lat/9411005 



Chiral Perturbation Theory and the 
Quenched Approximation of QCD* 

Maarten F.L. Golterman** 

Department of Physics 
Washington University 



> ; St. Louis, MO 63130, USA 

in 

o 
o 

^H 
^H 

ABSTRACT: The quenched approximation for QCD is, at present and in the foreseeable 
future, unavoidable in lattice calculations with realistic choices of the lattice spacing, vol- 
Q_i' ume and quark masses. In these lectures, I review the analytic study of the effects of 

quenching based on chiral perturbation theory. Quenched chiral perturbation theory leads 
to quantitative insight into the difference between quenched and unquenched QCD, and re- 
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veals clearly diseases which plague quenched QCD. A short review of the ideas underlying 
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chiral perturbation theory is included. 
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1. Introduction 

The lattice formulation of QCD has proven to be a powerful tool for computing QCD 
quantities of direct phenomenological interest, such as hadron masses, decay constants, 
weak matrix elements, the strong coupling constant, etc. (For reviews see for instance 
refs. [1,2], the proceedings of Lattice 93 [3], and the lectures of Rajan Gupta in these 
proceedings.) 

In order to perform such computations numerically, one obviously needs to consider a 
system with a finite number of degrees of freedom, which is accomplished by putting 
lattice QCD in a finite box. This box is then hopefully large enough to accomodate the 
physics one is interested in without serious finite size effects. This leads to the (minimal) 
requirement that the Compton wavelength of the particles of interest is sufficiently smaller 
than the linear dimension of the box, i.e. the mass has to be large enough for the particles 
to fit in the box. 

In order to have a small enough lattice spacing, small enough masses (in particular for 
the pion) and a large enough box size, one needs a large number of degrees of freedom in 
a numerical computation. It turns out that for QCD with realistic choices of the lattice 
spacing a, volume V and the quark masses (in particular the light quark masses), the 
presently available computational power is not adequate. The most severe problem comes 
from the fermion determinant, the logarithm of which is a very nonlocal part of the gluon 
effective action (specially for light quark masses). This nonlocality slows down currently 
available algorithms dramatically. 

In order to "circumvent" this problem, most numerical computations in lattice QCD have 
been done in the quenched approximation, in which one simply replaces the fermion de- 
terminant by one [4]. This amounts to ignoring all fermion loops which occur in QCD 
correlation functions (except those put in by hand through the choice of operators on the 
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external lines). While some handwaving arguments exist as to why this might not be un- 
reasonable, the quenched approximation does introduce an uncontrolled systematic error. 
Since the effect of a fermion loop is roughly inversely proportional to the fermion mass, this 
error is expected to be particularly large for quantities involving light quarks. Therefore 
one might expect that chiral perturbation theory (ChPT) is a useful tool for investigating 
the difference between quenched and unquenched ("full") QCD. 

In this talk, I will review a systematic approach to the study of the quenched approximation 
through ChPT [5,6,7,8]. There are two reasons why ChPT is useful in this context: 

• It turns out that ChPT can be systematically adapted to describe the low energy sector 
of quenched QCD [6] . It will therefore give us nontrivial, quantitative information on the 
difference between quenched and full QCD. 

• ChPT describes the approach to the chiral limit, and can be used for extrapolation of 
numerical results to small masses and large volumes. If these results come from quenched 
computations, one will of course need a quenched version of ChPT. (For finite volume 
ChPT, see refs. [9]. For quenched finite volume results, see refs. [6,7].) 

In this review, I will concentrate on the first point. First, I will give a short summary of the 
basic underpinnings of ChPT, in order to make these lectures more or less selfcontained. I 
will then give an early example of numerical results which can be understood using one loop 
ChPT. I will then go on to show how ChPT is developed for the quenched approximation, 
and use it for a quantitative comparison between full and quenched QCD. The quantities 
that I will discuss are fx/fn [6,8,10] and the octet baryon masses [11]. 

I will then address a number of theoretical problems that arise as a consequence of quench- 
ing. That such problems arise is no surprise, as quenching QCD mutilates the theory quite 
severely. It is however quite instructive to see what the actual consequences are. 
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2. A review of ChPT 

In this section I will summarize the basic ideas of ChPT, which is an efficient way to 
organize the information we can obtain from the QCD Ward identities for chiral symmetry. 
A partial list of useful references is [12,13,14]. Let us start with QCD with three flavors. 
The fermion part of the lagrangian is 

^QCD = QLP QL + QrP QR + QLMq R + q R M^q L . (1) 

q = (u, d, s) is a three flavor quark field, and the subscripts L and R denote the left- and 

righthanded projections qi R = ^(l ± 75) q. M is the quark mass matrix 

fm u \ 
M = m d . (2) 

V mj 

For M = 0, the lagrangian is invariant under the chiral symmetry group U(3)l x £/(3)^: 

q L ^U L q L , q R ^U R q R , U L ,U R eU(3). (3) 

Note that formally eq. (1) is invariant if we also let M transform as 

M -► U L MU R . (4) 

Axial U(l) transformations (i.e. U(l) transformations on the fields q^ and q R for which 
the phases are not equal) however are not a symmetry of the theory, since they are broken 
by the axial anomaly. The real symmetry group is therefore 5'L r (3)^ x SU(3) R x U(l) 
where the U(l) just corresponds to quark number conservation. It is furthermore believed 
that (still for M = 0) this symmetry is spontaneously broken to SU(3)y x U(l), where 
SU(3)y is the group of transformations for which Ui = U R G SU(3). The ensuing eight 
Goldstone bosons are the mesons ix , tt , K , K , K and rj, which transform in the octet 



representation of SU(3), denoted as the hermitian, traceless 3x3 matrix <fi: 

(5) 



f^ + JL n + K +\ 






These mesons acquire a nonzero mass due to the quark mass matrix M, as we will discuss 
below. 

One can now try to describe the physics of hadrons by constructing an effective lagrangian 
for composite operators with the quantum numbers of these hadrons. Explicitly, the 
operator H^j = qLiQRj (summed over color and spin indices, i and j are flavor indices) will 
couple to the Goldstone mesons, and one can introduce other operators for other hadrons, 
and write down an effective lagrangian of the form 

C eS (H,H\baryons,...). (6) 

This lagrangian, integrated over space-time, should be invariant under SU(3)l x SU(3)r x 
U(l) transformations H — > U^HUn etc. We can decompose H = RT, with R hermitian 
and positive, and E unitary. The spontaneous symmetry breakdown 577 ( 3) £ x SU(3)r x 
U(l) — ► SU(3)y x U(l) corresponds to R acquiring a nonzero vacuum expectation value 
rl. Dropping the other fields in £ e g-, the effective potential does not depend on E, and 
E = exp(2i(f)/f) can therefore be identified with the Goldstone bosons, eq. (5), leading to 
the nonlinear realization of chiral symmetry [15]. / is a constant with the dimension of a 
mass. (The overall phase of E corresponds to the rj which is not a Goldstone boson due 
to the anomaly [16,17], and has a mass of order 1 GeV.) The fluctuations of <fi describe 
the Goldstone mesons, whereas fluctuations in R around rl describe other heavy (scalar) 
mesons (which may or may not exist as narrow resonances in nature). 

At low energies (energies below the masses of any non-Goldstone hadrons) the effective 
lagrangian simplifies to 

£eff(£,M), (7) 

which has to be invariant under the transformations 

E - U L EUI 
M -► UlMU^. (8) 



This follows from the definition of H, which, ignoring fluctuations in R, is H = rE and 
from the invariance of the QCD path integral under the transformations eqs. (3, 4). Note 
that the unbroken symmetry SU(S)y is linearly realized (E — ► I/EI/i implies (j) ~* V(j)V '), 
whereas the broken symmetries are nonlinearly realized. 

Let us imagine that we can expand C e g in terms of derivatives and the quark mass matrix 
M. To second order in derivatives and linear order in M, (we will call this "0(p )") the 
most general expression is 

c2 



f 

C eS = J — tr(duEdutf) - v tr(M3$ + M+E) 



(9) 



where / and v are undetermined constants. Expanding C e Q to quadratic order in the 
meson field 4> (eq. (5)) using E = 1 -\ — j- + ... one will see that the kinetic terms have the 
standard normalization, and one can read off the tree level meson masses (for simplicity I 
will choose m u = rrid = m ) : 



2 8vm 2 Av{m + m s ) 2 9>v{m + 2m s 

' ~~ ~JT' m K = J2 ' m il = 3^2" 



m n = —FT' m K = To » m Ti - ^79 • ( 10 ) 



We see that the mesons acquire masses proportional to the square root of the quark masses, 
and in the chiral limit M — > they are massless as they should be. Note that for onshell 
mesons the combined expansion to second order in derivatives and to first order in M is 
consistent. Note also that we obtained our first nontrivial result: from eq. (10) it follows 

that 

2^2 1 2 /n\ 

m r) = z m K ~ 3 m 7n (11) 

which predicts a value of m^ about 3% too large. 

Let us proceed, and calculate scattering amplitudes A to some order in the loop expansion. 
In order to do this we need to introduce a cutoff A. The physical reason for this is that 
we are now ignoring all the hadronic physics at energies of order the p mass and beyond. 



We will not need to consider the Wess-Zumino term, see ref. [14] and refs. therein. 
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We therefore expect that we will have to choose A to be of the order of the p mass, i.e. 
or order 1 GeV . Using dimensional regularization and power counting one obtains for the 
contribution of a certain diagram to A [12,13] 

(2,)^(X»/V (j^f (})%( P 2 /A 2 ). (12) 

We identified A ks 4nf, which is of order 1 GeV (/ will turn out to be the pion decay 
constant, as we will see in a moment) [13]. Here pi stands for external momenta, and p 
denotes the square of any linear combination of these; it can also be a Goldstone meson 
mass squared. E is the number of external lines, and 

Ar = E^( rf - 2 )^+ L ' ( 13 ) 

d 
where L is the number of loops and V^ is the number of vertices with d derivatives (or 

d/2 powers of M, or any combination in between). F is a function containing logarithmic 

divergences in A arising from diagrams with loops. 

From eq. (12) we can draw the following conclusions. An amplitude A can be calculated as 
a series expansion in p /(4nf) if p and Wmeson are m uch smaller than A = (Airf) . In 
eq. (9) we only allowed 0(p ) terms, but we could have allowed higher order terms. From 
eq. (12) we see that treelevel contributions from 0(p ) terms come in at the same order as 
one loop terms from the 0(p ) lagrangian eq. (9): N — 1 corresponds to L — 1 and V4 = 
or L = and V4 = 1. Therefore the 0(p ) terms act as counterterms at the one loop 
level, and in fact need to be added at this order in order to absorb the cutoff dependence 
(cf. F in eq. (12)). A change in the choice of A can be absorbed by a shift in these 
counterterms. These conclusions generalize systematically to higher loops (L > 1). A very 
important conclusion of this analysis is that at any given order in p </(Anf) we only need 
a finite number of counterterms, and therefore the theory is predictive. Furthermore, the 
nonanalytic terms at L loops (contained in F) are determined by the 0(p ) lagrangian 
(i.e. by diagrams with Vj = for d > 2L). 
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As an example I will now discuss the decay constants f n and fx • First consider N = 0, so 
that the 0(p ) lagrangian at tree level is sufficient. The pion decay constant f n in QCD is 
defined by the matrix element of the lefthanded Noether current JLu, = ^Ll^ u L between 
a one pion state and the vacuum: 



i 



(0\j L ,(x)\7v + (p)) = -p fl f 7r e-^. (14) 

(It is this current which couples to the electroweak W bosons, through which the pion 
decays into a pair of leptons.) In ChPT this Noether current can be determined from eq. 
(9), leading to 

.•r2 f 

tr(T+9„SEt) = - J -duTT+ + ..., (15) 



JL„ = *-j- tr(T+a M EEt) = -ti,^ 



where T + is the appropriate 577(3) generator, and hence 

(0\jL^)\n + (p)}= l -p,fe-^. (16) 



We conclude that at this order f w — f (— 132 GeV) which determines the constant /. A 
similar calculation leads to 

fK = f = U, (17) 

which is a prediction, off by 22%. 

We see that the relation between f% and f^ is consistent with SU(3)y symmetry. The 
SU(3)y breaking introduced by the quark masses shows up at the one loop order (TV = 1). 
From ref. [14] 

f_K_ 1 5 mj /m|x _ 1 m\ /m|x _ 3 m\ ( m^, m\ - mj 

U + 416tt2/| S \A2) 2 16tt2/2 ° g I \2 J 4 16tt 2 /| V A 2 ' 16tt 2 /| 

(18) 

where L is a linear combination of coefficients of 0(p ) terms in the effective lagrangian 

[14]. With A = 1 GeV the experimental value fx/f-n — 1-22 can be reproduced with 

L = 1.1. Notice that the coefficient of log (1/A ) is proportional to m K — m\, so a change 



in the cutoff can be absorbed by a shift in L. We also see that the expansion parameters 
are 



m?T!~ om7 - Ifa% K009 ' (19) 

2 

which are small. In contrast, 1fi f. 2 rs 1.3, and hence the charm quark mass cannot be 
considered as small in the sense of chiral perturbation theory. 

An equation like eq. (18) can be used in two ways. First, we can determine L from the 
experimental values for fj{ and f n (fixing A at some value of order 1 GeV). This can then 
be used to predict other quantities, like /^ [14]. Second, eq. (18) can be used to fit results 
from numerical computations. Since these usually are performed at relatively large values 
of the quark mass (and hence the Goldstone meson masses), we can then employ results 
from ChPT to extrapolate to the physical values of quark masses. 

3. An example: B^ 

In this section, I would like to give an example of the definition and use of an electroweak 
operator in ChPT (for much more detail see e.g. refs. [13,18] and refs. therein). The kaon 
B parameter, B^ : which determines the strength of K — K mixing, is defined as 

\f K m\B K = (K°\O K \K ), (20) 

with 

Or = (sLl»d L )(s LllI d L ). (21) 

Ox transforms as a component of the (27, 1) representation of SU(3)i x SU(3)r. This 
symmetry property can be used to construct a corresponding operator in ChPT, which to 
0(p ) is uniquely given by 

0ChPT = l B f(d^) ds (d^) ds . (22) 



(Note that in both QCD and ChPT this operator is the product of two lefthanded currents.) 
The coefficient B is undetermined, and is therefore another free parameter in ChPT, like 
/ and v. To one loop order, we may calculate Bj{ from eqs. (20, 22), and we obtain (for 
degenerate quark masses, cf. ref. [19,7] for the nondegenerate case) 

1 - 6 £-s log -f + 0(p A ) contributions. (23) 

167T Z f4 A z J 

In ref. [7], besides Bjx, a different quantity By was studied, which, as we will see, has 
chiral one loop corrections larger than B^ has. For all details not discussed here I refer to 
ref. [7]. The quantity By is defined from 

\f\m\By = (K'°\[(s' a ^d' b )(s blf ,d a ) + (s , al „d' a )(s b ^d b )}\K°), (24) 

where a and b are color indices, d 1 and s' are new flavors of quarks, and K is a kaon 
built out of those. This is a technical trick to reduce the number of Wick contractions on 
eq. (24) to one. All quark masses are chosen equal. Note that I use a normalization for 
the meson decay constants which is different from that used in ref. [7]. With Bj± defined 
similarly with 7^ — ► 7^75 we have 

B K = By + B A . (25) 

One now can calculate By to one loop in ChPT, where in this special case (with degenerate 
quark masses) one can argue that the result is the same in the quenched and unquenched 
theories. The result is 

By = \B K - h\l 2 {m\) -h- \lh{m 2 K ) + 0(m 2 K ), (26) 

Z 4 m K o o 



him ) = -p / 77r-r A7 ^—-m - (27) 



where 

m f d^p 1 

J 1 J (2 7 r) 4 (p 2 + m 2 ) 1 

/3, 7 and 5 are coefficients of operators in ChPT which arise in the calculation of By to 
one loop (contrary to the case of Bj( there is more than one [7]). Bj± is given by the same 
expression with opposite signs for /?, 7 and 5. 
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Figure 1. Results for By versus m\ (= m K ) in lattice units (from, ref. [7]). 

From eq. (26) we see why it is interesting to consider By. Unlike in B^ there appears 
an "enhanced" chiral logarithm in By [20] - the second term on the righthand side of 
eq. (26). The extra factor m K in the denominator, if small enough, will enhance the size 
of this correction. If we repeat the calculation in a finite volume V = L , By becomes 
L dependent because the spatial part of the integral I2 now has to be replaced by a 
one loop momentum sum over momenta p = 2ivn/L with n E % (with periodic boundary 
conditions) . Volume dependence of this nature has actually been seen in numerical results. 
Fig. 1 is from ref. [7]. In this graph the points represent quenched numerical results. The 
solid lines show the result of a fit of the parameters ?>Bj( — |5, /3, and 7 (for a precise 
explanation of the fit see again ref. [7]). The graph shows that the numerical results are 
within errors consistent with one loop ChPT. Because the enhanced logarithms do not 
appear in B^ itself, one loop effects in this case are too small to be seen with the current 
statistical errors present in the numerical results for Bj{. 

This concludes our first example of a confrontation of one loop ChPT with numerical 
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computations. In this special case, the results from ChPT are unchanged by the effects 
of quenching. This is in general not true, and a systematic way of changing ChPT to 
correspond to the quenched approximation needs to be developed. I will do this in the 
next section, and then return to other, more recent, examples. 

4. Systematic ChPT for quenched QCD 

In this section I will outline the construction of a chiral effective action for the Goldstone 
boson sector of quenched QCD [6]. I will first introduce the formalism, and then show how 
it works in some examples. For early ideas on quenched ChPT, see ref. [21,5]. 

We will start from a lagrangian definition of euclidean quenched QCD. (We will restrict 
ourselves entirely to the euclidean theory which can be defined by a pathintegral. Hamil- 
tonian quenched QCD presumably does not exist.) To the usual QCD lagrangian with 
three flavors of quarks q a , a — u,d, s, we add three ghost quarks q a with exactly the same 
quantum numbers and masses m a , but with opposite, bosonic, statistics [21]: 

^quarks = ^Qaift + ™a)q a + ^q a (P + m a )q a , (28) 

a a 

where p is the covariant derivative coupling the quarks and ghost quarks to the gluon field. 

The gluon effective action produced by integrating over the quark- and ghost quarkfields 

vanishes, since the fermion determinant of the quark sector is exactly cancelled by that of 

the ghost sector. Note that the ghost quarks violate the spin-statistics theorem. Eq. (28) 

is the lagrangian definition of quenched QCD. 

We will now assume that mesons are formed as (ghost) quark - (ghost) antiquark pairs 
just like in ordinary QCD. This is basically equivalent to the notion that it is the dynamics 
of the gluons which leads to confinement and chiral symmetry breaking. The Goldstone 
particle spectrum of quenched QCD will then contain not only qq, but also qq, qq and qq 
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bound states. We will denote this 36-plet by 

*=(* *n~(^i q i)- (29) 

Note that the fields x an d X describe Goldstone fermions. 

The quenched QCD lagrangian eq. (28) with zero quark masses has a much larger sym- 
metry group than the usual U(3)l x U(3)r chiral group; it is invariant under the graded 
group U(3\3)i x U(3\3)r [6], where U(3\3) is a graded version of U(6) since it mixes the 
fermion and boson fields q and q. Writing an element U of £7(3 1 3) in block form as 

the 3x3 matrices A and B consist of commuting numbers, while the 3x3 matrices C and 
D consist of anticommuting numbers. 

We can now construct a low energy effective action for the Goldstone modes along the 
usual lines. We introduce the unitary field 

E = exp(2i$/f), (31) 

which transforms as E — > t/^EC/L with C/j, and IIr elements of £7(3 1 3). Because we are 
dealing here with a graded group, in order to build invariants, we need to use the supertrace 
str and the superdeterminant sdet instead of the normal trace and determinant, with [22] 

str(U) =tr(A) -tr(B), 
sdet(U) = exp(str log (£/)) = det(A - C B~ l D) / det(B) . (32) 

As one can easily verify, it is this definition of the supertrace that respects the cyclic 
property. To lowest order in the derivative expansion, and to lowest order in the quark 
masses, the chiral effective lagrangian consistent with our graded symmetry group is 

c2 



f 

C Q = J-str(dnZdn^) - v str(ME + MET), 



' \i^-"-> \i^ ) ~ u •"-' V-"^^ -r jviau j, (33) 
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where M. is the quark mass matrix 



m u 



M=[ 1 % m)> M= l ° m d ° I- (3-1) 

/ V m s 

As before, / and v are bare coupling constants which are not yet determined at this stage. 

The symmetry group is broken by the anomaly to the smaller group 
[SU(3\3)i j xSU(3\3)ji] (s)U(l) (the semidirect product arises as a consequence of the graded 
nature of the groups involved; the details are irrelevant for this talk). SU(3\3) consists 
of all elements U E U(3\3) with sdet(U) = 1. The anomalous field is $o — (V ~~ v')/V^, 
where the relative minus sign comes from the fact that in order to get a nonvanishing tri- 
angle diagram, one needs to choose opposite explicit signs for the quark and ghost quark 
loops, due to the different statistics of these fields, rj is the field describing the normal r/ 
particle, while fj' is the ghost rj consisting of ghost quarks and ghost antiquarks. We will 
call the field $q the super-?/ field. The field $q °^ s ^ r l°g ^ = l°g s det E is invariant under 
the smaller symmetry group, and we should include arbitrary functions of this field in our 
effective lagrangian. In analogy to ref. [14], the correct chiral effective lagrangian is 

C =V 1 (^ )str(d^d fl ^) - V 2 ($o)str(ME + MT^) 

+ Vo(<S>o)+V 5 (<S>o)(d^ ) 2 , (35) 

where the function multiplying i str (A4 E — .ME * ) can be chosen equal to zero after a field 
redefinition. This lagrangian describes quenched ChPT systematically, as we will show 
now with a few examples. 

For our first example, let us isolate just the quadratic terms for the fields rj and fj' ', choosing 
degenerate quark masses for simplicity. We expand 

Vi(*o) = y + ---, 

y 2 ($ ) = v + ..., 

V ($ ) = constant + /j 2 $q + ---, (36) 

y 5 ($ ) = « + ..., 
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and obtain 

Ctf, ff) =\(d^') 2 - l -(d^>f + i«( V - drff 

+ ^rnliv') 2 - \ml{fif + ±n\rf - fj'f + ..., (37) 

where m 2 = 8mv/f. The relative minus signs between the rj and ff terms in eq. (37) 
come from the supertraces in eq. (35), and are related to the graded nature of the chiral 
symmetry group of quenched QCD. 

The inverse propagator in momentum space, 

(p 2 + ml)(l _° 1 )+(^ 2 + «p 2 )(_ 1 1 - 1 ), (38) 

clearly cannot be diagonalized in a p independent way, which is quite different from what 
one would expect from a normal field theory! Treating the \i + ap term as a twopoint 
vertex, one can easily show that the repetition of this vertex on one meson line vanishes, due 
to the fact that the propagator matrix I „ _i I multiplied on both sides by the vertex 

( 1 -l\ 

matrix I _- -.1 gives zero. This result coincides exactly with what one would expect 

from the quark flow picture for rj propagation, as depicted in fig. 2. The straight-through 
and double hairpin contributions do not contain any virtual quark loops, and are therefore 
present in the quenched approximation. All other contributions should vanish because 
they do contain virtual quark loops, and this is exactly what happens as a consequence 
of the (admittedly strange) Feynman rules for the propagator in the 77-77 sector! This 
propagator is given by the inverse of eq. (38) and reads 

1 fl \ fi 2 + ap 2 (\ i\ 



p 2 + m 2 \0 — 1 J (p 2 4- m 2 ) 2 \1 1 

in which the two terms correspond to the two first diagrams in fig. 2. 
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Figure 2. Quark flow diagrams for the r( propagator in full QCD. 

From eq. (39) we learn several things. First, because \x , which in full ChPT would 
correspond to the singlet part of the rj mass, appears in the numerator, we need to keep 
the rf (and its ghost partner) in quenched ChPT: it cannot be decoupled by taking /j, 
large. Second, this "propagator" is definitely sick, due to the double pole term. It should be 
stressed here that this double pole term is an unescapable consequence of quenched QCD, 
and does not result from our way of setting up chiral perturbation theory. In the case of 
nondegenerate quark masses, this double pole also shows up in the it and r\ propagators, 
due to mixing with the rj . I will return to these strange properties of quenched QCD in 
section 6. 



Figure 3. One loop pion selfenergy in quenched ChPT. 

As a second example, we will consider the (charged) pion selfenergy at one loop, again with 
degenerate quark masses. I will set a = for simplicity. At one loop, the pion selfenergy 
only contains tadpoles, with either (j) or x lines (cf. eq. (29)) on the loop. Also, on the <p 
loop, one can have an arbitrary number of insertions of the vertex /j, if the internal <f> line 
is an SU(3) singlet. These various contributions are drawn in fig. 3, where a solid line 
denotes a <j> line, a dashed line denotes a x line, an d a cross denotes a /j, vertex. One finds 
that the diagrams with the <p and x lines on the loop without any crosses cancel, and then, 

16 



of course, that the diagrams with more than one cross vanish, using our earlier result for 
the r)'-fj' propagator. We are left with only one term, and the result is 

_ . , 2ml (J 2 f d 4 k 1 . t . 

*m = 1 ?tJt2^wt^- (40) 

The pion selfenergy is logarithmically divergent, but this nonanalytic term is completely 
different from those that arise in the unquenched theory, as it is proportional to fj, . One 
can easily convince oneself that the diagrams in fig. 3 which cancel or vanish correspond 
to diagrams with virtual quark or ghost quark loops in the quark flow picture. (For early 
discussions of the quenched pion selfenergy in the quark flow picture, see refs. [21,5].) 

Before I go on to look at some quantitative results, I would like to discuss one aspect of 
the chiral expansion in quenched ChPT. The chiral expansion is basically an expansion 
in the pion mass (see e.g. ref. [12]). However, as we have argued above, in quenched 
ChPT there is unavoidably another mass scale, namely the singlet part of the rf mass, 
fj, . For our expansion to be systematic as an expansion in the pion mass, we would 
have to sum up all orders in \i , at a fixed order in the pion mass. This is clearly a 
formidable task. In order to avoid this complication in a systematic way, we can think of 
fj, /3 (which turns out to be the natural parameter as it appears in the chiral expansion) 
as an independent small parameter. To check whether this makes any sense, one may note 
that from the experimental value of the 77' mass one obtains a value \i /3 ~ (500 MeV) , 
which is roughly equal to the kaon mass squared, m K . Of course, for quenched QCD 
the parameter fj, need not have the same value, after all quenched QCD is a different 
theory. A lattice computation of this parameter [23] gives Mouenched/^full ~ 0.75. (a can 
be estimated from 77-77' mixing, and is very small.) Finally, one may also note that both 
fj, and a are of order 1/N C , where N c is the number of colors [17]. I will return to this 
point in section 6. 
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5. Quantitative comparison of quenched and full ChPT 

Let us first consider the quenched result for the ratio of the kaon and pion decay constants 
fx and f n [6,8]. I will set a = and take m u = m^ = m: 



Jk J quenched l" 7r Tn 



™i' l 0g ( 2 4L _ ! | _ , 



2(m K - ml) & I raf 



2 2 

+ -t?]f L - (41) 



'7T 



L is again a certain combination of 0(p ) couplings of the quenched chiral lagrangian, which 
does not have to be equal to L in eq. (18), since it comes from a different theory (quenched 
QCD). (For example, L is cutoff dependent, while L is not.) Because of this, the result 
eq. (41) is not directly comparable to the equivalent result for the full theory. In other 
words, in order to compare quenched and full QCD, we have to consider quantities which 
are independent of the bare parameters of the effective action. (Alternatively, we would 
need to extract the values of the bare parameters from some independent measurement or 

lattice computation, in this case, we would need independent determinations of L in both 

1/3 4/3 
quenched and full QCD.) In the full theory, (f^f^ )/fi{ * s such a quantity [14], but in 

the quenched theory this quantity is not well defined, due to the double poles which occur 

in the propagators of neutral mesons. 

We will therefore choose to consider a slightly different theory, in which sufficiently many 
charged (i.e. off-diagonal) mesons are present [8]. This theory is a theory with two light 
quarks m u = m^ = m and two heavy quarks m s = m s > = ml ' . This theory contains a ud 
pion 7T, an ss pion n and a us kaon K, with (tree level) mass relation 

™%= ^™l + ™l>)- (42) 



One can show that the ratio fK/V7nfn' ^ s independent of the low energy constants L. 
For the quenched theory we find 



. 1-loop , 



\[~hU> J minTnrhpH 16yr f* 



quenched 
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ml + m "' logf^l-1 



2 



2(mp-ml) ' \ml 



(43) 



whereas in the full theory 

,2 \ / Im 2 



, 1— loop 
f'K \ 1 1 



2 / m K I 2 / m K 

m * l0g | ~^J ) + "V l0g | _ T 



v?^; full " * *&% "■* "* ■ =*■ " '"*' '" s ' ^ ' ' (44) 

Note again that the logarithms in the quenched and unquenched expressions are completely 
different in origin. 

We may now compare these two expressions using "real world" data, where we'll determine 

the value of the n mass from the mass relation eq. (42). With m n = 140 MeV, m^ = 

494 MeV and y 2 /3 = 0.75 x (500 MeV) 2 we find 

\ 1— loop 

Ik \ 
h ' 1.049, 



V JKJn / q UCnc hed 
\ 1-loop 
4== = 1-023, (45) 

V/TT-Wfull 

a difference of 3%. If we choose /j, /3 = (500 MeV) , we find a difference of about 4%. 
This difference is small. Note however, that this is due to the fact that for this particular 
ratio, ChPT seems to work very well, both for the full and the quenched theories. If one 
only considers the size of the one loop corrections (the numbers behind the decimal point), 
the quenched and full results are very different. It is also possible, and in fact not unlikely, 
that part of the difference between the full and quenched theory gets "washed out" by the 
fact that we are considering a "ratio of ratios" . It follows that the relative difference is a 
lower bound on the difference between the quenched and full values of the decay constants. 
For another quantity for which the difference between quenched and full ChPT has been 
calculated, see ref. [8]. 

Recently, numerical results for quenched fx/fn have become available which are precise 
enough to make a comparison with eq. (43) interesting. These results are shown in fig. 4 
[10], in which the solid line depicts a fit of eq. (43) to the numerical results. The quantity 

X is defined as 

mi + mi/ i ml, , 

x = ^r JL 9 =07 log -7- - 1 . ( 46 ) 



2(mp - m 2 ,) \ m n 
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Figure 4. fj^/ifirfn') versus X (see text for explanation) (adapted from ref. [10]). 



and what was fitted is 



with 



[K 

firU' 



1 + SX, 



^ 3 /3 



(47) 



(48) 



The solid line corresponds to 5 = 0.133, and the dashed lines correspond to a variation of 
±0.033 in this value for 6. We expect 6 = 0.137—0.182 from eq. (48) for /j?/S = (0.75- 
1) x (500 GeV) . If one plots the same quantity using eq. (44) instead of eq. (43), the 
results look much less convincing, however, the errors in the numerical results are still 
quite large. At this stage one can conclude that numerical results are consistent with one 
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loop quenched ChPT. It would be interesting to have numerical results with smaller errors. 
For more detail, see ref. [10]. 

Next, I will review some work on baryons in quenched ChPT by Labrenz and Sharpe [11]. 
They calculated the one loop corrections to the octet baryon mass coming from the cloud 
of Goldstone mesons. They employed an effective lagrangian for quenched heavy baryon 
ChPT, constructed using the same techniques as described in section 4. In the case of 
degenerate quark masses, the result for the nucleon mass is 



m^v —m — tt(D — 3F) c. m-n + 2(6p — ?>bp)m 1T 



+ 



2 V ' SrrVf 

|(D - 3F)(2£> + 3 7 ) + ^a(D - 3F) 2 



(49) 



In this equation, m, D, F, bp>, bp and 7 are bare parameters which occur in the baryon 
effective action, m is the bare "average" octet mass, D and F are the well known baryon- 
meson couplings, bp, and bp are low energy constants which arise as a consequence of 
renormalization (see for instance refs. [24,25]). 7 is a new coupling which occurs because 
of the unavoidable presence of the super-?/ in the quenched approximation. The term 
proportional to fj, comes from a diagram with a cross on the (f) internal line, i.e. an insertion 
of the n twopoint vertex. Note that in this case there are also one loop corrections not 
involving \x which survive the quenched approximation, in contrast to the pion selfenergy, 
eq. (40), or fx/firj e Q- (41). The authors of ref. [11] then calculated the coefficients using 
full QCD values for the various parameters (from ref. [25]). With a = and 7 = (7 = 
is consistent with available information, which however is limited [26]), they obtained 

m N = 0.97 - 0.5 — m^ + 3.4mJ - 1.5mJ, (50) 

with 6 as in eq. (48) and 5 ~ 0.2 for the full theory. 

In ref. [11], eq. (49) was also compared to recent numerical results from ref. [27]. These 
data are presented in fig. 5, where the scale a = 1.63 GeV is set by f n [11]. If one 
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calculates the coefficients in eq. (49) by "fitting" the four data points, one finds 

m N = 0.96 - I.Otott + 3.6ml ~ 2.0mJ. (51) 

This is only four data points for four parameters, and the "fit" is very sensitive to for 
instance an additional m\ term. Nevertheless, from the agreement between eq. (50) and 
eq. (51) it appears that it is reasonable to apply ChPT to the results of ref. [27]. Note 
that the individual terms in eq. (50) are quite large for the two higher pion masses in 
fig. 5 (this is not unlike the case of unquenched ChPT). From fig. 5 it is also clear that 
(m/fir) quenched ¥" (jn/ fn) full because of the term linear in eq. (49), which is absent in full 
ChPT. 

Labrenz and Sharpe then went on to consider octet mass splittings. In order to remove 
effects which can be accomodated by a change of scale, they calculated the ratios 

Til ' 

Rij = — , i,j = N,A,V,E (52) 

■' rrij 

in quenched ChPT, and compared these with similar ratios obtained from ref. [25]. They 
assumed that all bare parameters in the equations for the octet masses (for explicit ex- 
pressions, see their paper) are equal in the full and quenched theory, and then calculated 
the ratios 

D quenched 

Ti J = Dfull • ( 53 ) 

R ij 

With the assumption that the bare parameters of the quenched and full theories are equal, 
bp and bp drop out of the ratios, and with 7 = 0, a = and D and F equal to their full 
QCD values, they obtain 

r^ N = 1 + 0.19(5/0.2) + 0.13 = 1.31(1.27] for 5 = 0.2[0.15], 

r~N = 1 - 0.46(5/0.2) + 0.43 = 0.97(1.09] for 8 = 0.2(0.15], (54) 

r AN = 1 - 0.39(5/0.2) + 0.26 = 0.87[0.97] for 5 = 0.2(0.15]. 

(The choice 5 = 0.15 corresponds roughly to the value reported in ref. [23].) 
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Figure 5. TTie nucleon mass from the lattice [27] (copied from ref. [11]). The curve is 
from a fit to the form m^ =m + am^ + brn^ + cm^ . 

From this, one would conclude that one can expect errors from quenching of at least 20% 
in the octet splittings. These differences between the quenched and full theories cannot be 
compensated for by a change in scale between quenched and full QCD. 

At this point I would like to comment on the above mentioned assumption that was used 
in order to obtain eq. (54). Let us consider in particular the parameters bjj and bp. They 
correspond to higher derivative terms in the baryon-meson effective action, and are needed 
in order to absorb the UV divergences which arise at one loop in ChPT. Since the size 
of these divergences is in principle different between the full and quenched theories, one 
expects that ^D^uenched an d ^F, quenched can be different from &£} f u n and bjrf u n. If we 
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want to proceed without assuming that the quenched and full 6's are equal, we have to 
consider ratios of quantities independent of the parameters bp, and bp. The situation is 
essentially the same as in the case of fx/ In- From the available results [11], only one ratio 
independent of bp and bp can be formed: 

y = ^4_ (55) 

If we expand Y in the Goldstone meson masses using ChPT, Y — 1 measures the deviation 
from the Gell-Mann-Okubo formula (cf. ref. [25] for the full theory). 

Setting m^ = keeping only m^ as in ref. [11], one finds 

3 
Quenched = 1 " 1T046^|^ (d 2 - 3F 2 ) + 1.33335^ (d 2 - 3F 2 



in 



3 



Y Ml = 1 - 0.4125^-^2 [D 2 - 3F 2 ) . (56) 

(The parameter 7 drops out of this particular combination and we again take a = 0.) 
These quantities still depend on the other bare parameters, D, F and m. Again, they 
could be different in the quenched and full theories, and I will leave the quenched values as 
free parameters. Substituting mp; = 495 MeV, f n = 132 MeV, -Df u ii = 0.75, Ff u ^ = 0.5 
and mf u n = 1 GeV [25] finally gives 



Vcnchcd =1 _ 00214 + 



Y t 



full 



5 

0.293 0.306 

0.2 



/ 7~)2 OTp2\ 

\ 1 quenched 

^quenched/-'- ^- Te > / 



(57) 



For any reasonable values of m and 5, and for (D — 2>F ) , , not too far from its full 

theory value of —0.1875, the difference between the quenched and full theories as measured 
by the ratio ^quenched /^full ^ s n °t more than a few percent. Of course the same comment 
that applied in the case of fx/ fir applies here, that part of the difference may have been 
washed out by taking "ratios of ratios" . Summarizing, the conclusion of this analysis seems 
to be that the error from quenching for octet baryon masses is at least a few percent, and 
could be as much as 20%. 
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6. A sickness of quenched QCD 

Let us again consider the quenched result for fx/f-K-, ec l- (41), as a function of the quark 
masses (using treelevel relations between meson masses and quark masses), 
1-loop l( 3 

quenched 



j^ / aucnehed ^ n jit 



log 1 



+ L — term. 



_2(m s - m u ) m u 

From this expression it is clear that we cannot take m u — » keeping m s fixed, or, to 
put it differently, that if we take both m u and m s to zero keeping the ratio fixed, the 
limit depends on this ratio, and is not equal to one! This is quite unlike the case of full 
ChPT, where one can take any quark mass to zero uniformly, and deviations from SU(3) 
symmetry due to this quark mass vanish in this limit. Technically, the reason for this 
strange behavior is that there is another mass fj,, which, as we argued before, cannot be 
avoided in quenched ChPT. This mass is related to the singlet part of the rj mass, and is 
not a free parameter of (quenched) QCD. Even if we do not consider any Green's functions 
with rj external lines, this mass shows up through the double pole term in eq. (39) on 
internal lines. Because of the double pole, such contributions can lead to new infrared 
divergences in the m^ — ► limit. This problem with the chiral limit of quenched ChPT 
shows up in other quantities, such as meson masses and {ipip) [6,7]. 

A question one might ask is whether this problem is an artifact of one loop quenched 
ChPT [8]. For instance, if we would sum all contributions to the r\ propagator, maybe the 
double pole term would become softer in the p — ■> limit, improving the infrared behavior 
of diagrams in which the double pole terms appear. Let us address this question in the 
chiral limit, m a = 0, where the problem is most severe. In the full theory, we can write 
the fully dressed rj' propagator as 

Z(P » (58) 



p 2 + S(p)' 
and define Hp(p) = S(p), which onshell is the rj mass in the chiral limit. Likewise, in the 
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quenched theory we can write the dressed rf , r)' propagator as 

Zq( P ) ' ' u ' '' " ' ' 



(59) 



p 2 {O -1) (p2)2 \l 1 

which defines Hn(p)- To leading order in l/iV c , these two definitions of /U (p) should lead 
to the same result: 

^(p) = m|-(p)(i + o(^)). 

We also believe that jJbj?{p = 0) is not equal to zero, since we expect the rf in the full 
theory to remain a well-behaved, massive particle in the chiral limit. This implies, in sofar 
as we can rely on the large N c expansion, that Hq(0) ^ 0, and that the double pole in 
eq. (39) is a true feature of the theory. (The argument can be repeated at nonzero quark 
masses, which is necessary if the chiral limit of quenched QCD does not exist.) 

While this argument is not very rigorous, I believe that the foregoing discussion implies 
that the chiral limit of quenched QCD really does not exist. This belief is furthermore 
supported by the following remarks: 

• Sharpe [7] has summed a class of diagrams in the case of degenerate quark masses for a 
very simple quantity (the pion mass), and found a result that is actually more divergent 
than the one loop result. 

• With nondegenerate quark masses there are many more diagrams that are infrared 
divergent in the chiral limit, and it is even less probable that resummation will improve 
the situation. 

• Any mechanism improving the infrared behavior would have to work for each divergent 
quantity. One expects that such a mechanism would be related to the double pole term in 
the rj propagator, which created the problem in the first place. But this seems unlikely in 
view of the arguments given above. 

• The bare quark mass parameter appearing in the chiral effective action is not the same 
as that appearing in the (unrenormalized) QCD lagrangian. But one can argue that the 
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two bare quark masses should be analytically related, and the infrared problem is not just 
a problem of quenched ChPT, but of quenched QCD. 

7. Conclusion 

The quenched approximation leads to an unknown systematic error in all lattice calcula- 
tions that employ this approximation. It would of course be very nice to have a parameter 
that interpolates between full and quenched QCD, and in principle the quark masses could 
play such a role, since one expects that quenched QCD corresponds to full QCD with 
very heavy quarks. One would have to distinguish here between valence and sea quark 
masses, and it is the sea quark mass that would play the role of such a parameter. This 
distinction can indeed be made by considering so-called partially quenched theories [28], 
but no practical scheme to implement this idea is known. 

Quenched QCD can be defined from a euclidean pathintegral as rigorously as full QCD. 
In this talk I have explained that euclidean quenched ChPT can be used as a tool for a 
systematic investigation of quenched QCD. Quenched ChPT does not quite accomplish a 
task equivalent to that of an interpolating parameter. Since the bare parameters appearing 
in the quenched and full chiral effective actions are not the same (as explained in section 5) 
one cannot directly compare quantities calculated in full and in quenched ChPT. However, 
one can calculate combinations of physical quantities which do not depend on the bare 
parameters, and in that case a direct comparison between quenched and full QCD is 
possible, as we demonstrated with an example involving meson decay constants. This 
makes it possible to estimate lower bounds on the differences which come from quenching; 
these estimates are dependent on the values of the meson masses, which can be taken to 
be the (known) independent parameters of the theory. For realistic values of these masses, 
such differences turn out to be of the order of a few percent for ratios of decay constants 
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and for baryon octet splittings. 

The disadvantage of this more conservative approach is that part of the difference maybe 
hidden, because these specific combinations of physical quantities maybe less sensitive to 
the effects of quenching than other quantities of interest. This is particularly clear in 
the example of baryon octet masses. In this case, a comparison based on the assumption 
that the bare parameters of the full and quenched effective theories are the same, lead 
to differences of up to 20% and more. Of course, it is not known to what extend this 
assumption is valid. 

The differences between the quenched and full theories become markedly larger for decreas- 
ing quark masses. This is due to the fact that new infrared divergences occur in quenched 
QCD, which do not have a counterpart in full QCD. These divergences lead to the nonex- 
istence of the chiral limit for quenched ChPT (as discussed in section 6). The origin of this 
phenomenon can be traced to the special role of the rj in the quenched approximation. 
In the quenched approximation, the rj is a Goldstone boson (it develops massless poles in 
the chiral limit), but an additional double pole term arises in its propagator, rendering it 
a "sick" particle. For nondegenerate quark masses this problem is also inherited by the n 
and the r\. In section 6 I argued that the nonexistence of the chiral limit is a fundamental 
feature of quenched QCD. 

In principle therefore, the chiral expansion breaks down for quenched QCD. For very small 
quark masses, at fixed values of the singlet part of the rj mass fj, , the expansion becomes 
unreliable. In order to make progress, one may take the expansion to be an expansion in 
fi /3 (which was shown to be roughly equal to m K phenomenologically), with coefficients 
which are functions of the quark mass. These functions sometimes show divergent behavior 
in the chiral limit {e.g. the one loop correction to fx/ fir)- If such divergences occur, the 
expansion is only valid for a range of quark masses which are neither too small, nor too 
large. It would be interesting to see whether this point of view can be made solid. One 
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way to check whether quenched ChPT makes any sense, is to compare its predictions with 
numerical results. Results for fx/fn (section 5) and By (section 3) are consistent with 
one loop quenched ChPT. 

It is in principle interesting to study any quantity which is being computed in quenched 
lattice QCD in ChPT, for those quantities for which ChPT is applicable (meson masses, 
decay constants, condensates and the kaon B parameter have been calculated [6,5,8]). As 
discussed, this includes not only Goldstone meson physics per se, but also chiral corrections 
to baryon masses [11], and for the same reason, to mesons containing heavy quarks. 

Recently, also attempts have been made to compute pion and nucleon scattering lengths 
[29,30] from quenched lattice QCD. If one tries to calculate the 1 = pion scattering 
amplitude in quenched ChPT, one actually finds that the imaginary part is divergent at 
threshold, even for nonvanishing pion mass [31]! Again, this can be related to double 
pole terms in the rj propagator. Apparently euclidean quenched correlation functions in 
general cannot be analytically continued to Minkowski space-time. (The euclidean four 
pion correlation functions are well defined.) Further work is needed on pion scattering 
lengths. 
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